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1. Introduction 

Particles hopping on a one-dimensional lattice with hard-core exclusion interactions 
provide a simple framework for the study of interacting many-body systems [HE]- In 
particular, when the particle hopping is asymmetric, a macroscopic particle current 
results and the system attains a nonequilibrium steady state (NESS) in which detailed 
balance is not satisfied (see [3] for a recent review). 

A fundamental example of such a system is the Asymmetric Simple Exclusion 
Process (ASEP) [H [5]. Here particles attempt hops to the right neighbour site with 
unit rate and to the left neighbour site with rate q. The hop is carried out when the 
destination site is empty. The special case q = is known as the Totally Asymmetric 
Simple Exclusion Process (TASEP). With periodic boundary conditions the NESS of 
the ASEP has a very simple form (all allowed configurations of particles are equally 
likely) yet dynamical properties such as diffusion of a tagged particle P, [2, El E], 
and large deviations of the current [10] have proved to be non-trivial. When open 
boundary conditions, where particles attempt to enter and exit at the left and the right 
boundary, are used instead of periodic, the NESS takes on a non-trivial form. It may be 
represented by a matrix product state [H] in which the steady state probabilities for each 
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configuration are obtained from products of two matrices D and E according to whether 
each site is occupied or empty in the configuration. These matrices obey a quadratic 
algebra which provides a motif from which all the steady state probabilities may be 
generated. This quadratic algebra is related to the q-deformed harmonic oscillator 

nana ei us HE]. 

A generalization of the ASEP is to the case of two species of particles. A well 
studied model is that of usual (first-class) and second class particles. In this case both 
first and second-class particles hop to the right with unit rate and to the left with rate 
q. However if the site to the right of a first-class particle is occupied by a second-class 
particle the first and second-class particle exchange places with rate 1. Conversely, if 
the site to the left of a first-class particle is occupied by a second-class particle they 
can exchange with rate q. Thus a second-class particle behaves as a hole from the 
point of view of a first-class particle but behaves as a particle from the point of view 
of a hole. When q is lower than 1 the second class particles will move forwards in an 
environment of a low density of first class particles but backwards in a high density 
environment, therefore the introduction of a second class particle is a useful tool in the 
study of the microscopic structure of shocks p21 [181 EE EDI EE) • The stationary state of 
a periodic system containing second and first-class particles has been obtained using the 
matrix product formulation by Derrida et al. |20j. Here the three matrices D, E and A 
(corresponding to first-class particles, second-class particles and holes respectively) obey 
a quadratic algebra closely related to that of [11]. Quadratic algebras for two species 
exclusion process have been further explored [221 E31 El] and a general classification 
of quadratic algebras has been made in [25]. The range of matrix product states and 
quadratic algebras relevant to exclusion process and other stochastic systems have been 
summarized in [3]. 

A natural generalization of the two species case of first and second class particles 
is to the multispecies process where there is a hierarchy amongst the different species. 
That is, for a system with N species, the iVth class particles are treated by all other 
classes of particle as holes, (iV— l)th class particles treat iVth class particles as holes but 
are themselves treated as holes by first class, second class,.... (N— 2)th class particles, 
and so on, up to the first class particles which treat all other species as holes. We refer 
to this model as the iV-ASEP. In the physics literature, the N = 3 totally asymmetric 
case was considered by Mallick, Mallick and Rajewsky [26] and a matrix product steady 
state was determined. It was shown that the matrices obeyed more complicated relations 
than the quadratic algebra of the two species case. 

In the probabilistic literature Ferrari and Martin [27] provided a construction for 
the A-TASEP whereby the dynamics is related to the dynamics of an N line process 
i.e. N coupled TASEPs. The steady state for the N line process has each of its 
configurations equally likely. Therefore, to sample the configurations of the A-TASEP 
according to their steady state probabilities one picks from a uniform distribution an 
N line configuration and projects this onto the corresponding iV-TASEP configuration. 
Ferrari and Martin also showed that the iV-line construction could be interpreted as a 
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system of N queues in tandem with priority customers. 

In a recent work [28J it was shown how the construction of Ferrari and Martin can 
be inverted and the steady state probabilities written down as matrix products. A proof 
was given for the case N = 3 and the matrices for general N were written down in a 
hierarchical fashion. 

In the present work we generalize the solution to the partially asymmetric case and 
provide the full proof of the matrix product state for arbitrary N. Let us summarize 
the key points of our solution. The matrices Xj , appearing in the matrix product 
expression for the steady state (1241) . depend on J, the species present at the site to 
which the matrix corresponds, and N the number of species represented in the system. 
The matrices are defined in eq. (1251) in a hierarchical way: the matrices for an iV 
species system are expressed in terms of those for an N— 1 species system. This in 
turn means that the weights of the iV species system may be expressed in terms of the 
weights of an N— 1 species system via a transfer matrix which is defined in eq. (144)) . 
For the two species case (N = 2) the matrices obey a quadratic algebra ( [1011111)121) 
and this allows reduction relations which relate the weights of a system of L sites to 
those of a system of L — 1 sites. However, for N > 3 the quadratic algebra is replaced 
by a more complicated set of relations (I69ll70|l7ip involving additional 'hat' matrices, 
again defined hierarchically in equations (I73)74p . Both the matrices Xj N ^ and the hat 
matrices Xj^ are expressed recursively in terms of auxiliary matrices defined in 



equations (T27H33]) . These auxiliary matrices are themselves tensor products of the four 
fundamental matrices 1, 5, e, A which appear in the N = 2 solution. The algebraic 
properties of the auxiliary matrices, in particular the symmetry relation (15?]) and the 
commutation relations ( I6TH631 . are the key to the proof of the quadratic relations obeyed 



by the matrices Xj N ^ and Xj N ^ which, in turn, furnish the proof of the matrix product 
representation of the stationary state. 

The transfer matrix mentioned above allows us to investigate how the construction 
of Ferrari and Martin generalizes to the partially asymmetric case. For TASEP, the 
transfer matrix implements the Ferrari-Martin construction explicitly. However, in 
the system with partial asymmetry, the queueing process interpretation does not hold 
anymore and is replaced by a more general recurrence between systems with N and 
N — 1 species. 

The paper is structured as follows: in section [21 we define the multispecies 
asymmetric exclusion process. In section [3], we write the matrix product representation 
of the stationary state of this model. In section H] we give an interpretation of this matrix 
product representation in terms of a transfer matrix. Then, after writing the algebraic 
relations obeyed by the auxiliary matrices related to the matrix product representation 
in section [3J we give a complete proof of the matrix product expression in section 



El |Appendix A| is devoted to the calculation of some traces of product of matrices, 



Appendix B proves a classification of nonzero elements of the transfer matrix, while 



Appendix C contains the proof of a special case of the identities required to ensure the 



matrix product expression is valid. 
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2. Multispecies ASEP 

We consider the multispecies asymmetric simple exclusion process with both forward 
and backward jumps on a one dimensional lattice with periodic boundary conditions. 
This stochastic model is defined on a configuration space such that each of the L sites of 
the lattice can be occupied by at most one particle (exclusion rule). Each particle has a 
label which is an integer between 1 to N, the 'class' of the particle. (We use the terms 
'class' and 'species' interchangeably.) The unoccupied sites (holes) will be considered 
as particles of class 0. The stochastic dynamics can be expressed in terms of exchanges 
of particles at neighbouring sites. The transitions which can occur depending on the 
classes of both particles are 

JK -> KJ with rate 1 if 1 < J < K < N (1) 

KJ -> JK with rate q if 1 < J < K < N (2) 

JO -> 0J with rate 1 if 1 < J < A (3) 

J -> JO with rate q if 1 < J < A . (4) 

All classes of particles jump to the right with rate 1 and to the left with rate q if the 
destination site is empty. When two particles of different class are on neighbouring sites, 
they can exchange with rate 1 if the particle with the smallest class is on the left, and 
with rate q if it is on the right. Thus, for a particle of class r, all the particles of class 
larger than r behave as holes. At this point, it might seem natural to consider holes as 
particles of class A + 1 rather than 0. However, we do not adopt this convention as it 
would make the expression of the stationary state more complicated in the following. 

We use the site variable r% — 0, 1, . . . , A. If T\ = the site is empty; if T/ = r > 
site / contains a rth-class particle. Let us denote by r = (n, . . . , Tl), a configuration of 
the system. The dynamics of the system can be encoded in a Markov matrix M. The 
time evolution of the probability Pt(r) to be in a configuration f at time t is given by 
the master equation 

|p t (r) = ^Af(f,f')P t (f'). (5) 

f' 

The matrix M is a (N + 1) L by (N + 1) L matrix which acts on the configuration space. 
As the numbers of particles of each class are conserved by the dynamics, we will restrict 
ourselves to a configuration space with fixed number of particles and holes. We call P r 
the number of particles of class r. The restricted configuration space Q(P Q ,Pi, ...,P^) 
has dimension 

and the restricted Markov matrix which acts on it is by 
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3. Matrix product formulation of the stationary state of the iV-species 
ASEP 

The matrix product formulation was first used to solve the TASEP on a lattice of 
length L with open boundary conditions [IT] . It was then extended to the 2-ASEP on 
the ring Z^ [20], which is our starting point. In this case the site variable is T\ = 0, 1, 2 
which implies that site / is empty, contains a first-class particle or contains a second- 
class particle, respectively. For the 2-ASEP the matrix product formulation consists 
of writing the weight of a given configuration as a trace of a product of L matrices 
corresponding to the classes of particles on the different sites. In the matrix product 
formulation [20] it has been proved that the stationary measure may be written as 

1 



■W(t) 



y ' Z{P ,P 1 ,P 2 ) 
where the weight of the configuration is given by 

L 



W{f) = Tr 



i=i 



(7) 



(8) 



and Tr means the trace of the product of matrices X n . The normalization Z(Pq, Pi, P2) 
is chosen so that the sum of all the probabilities is equal to 1 i.e. Z(P , Pi, P 2 ) is the 
sum of weights for configurations with the correct numbers Pq, Pi, P 2 of holes and each 
species of particles. The matrices X Tl are given by 

X = E, X 1 = D, X 2 = A, (9) 

that is: if the site is empty we write a matrix E; if the site contains a first-class particle 
we write a matrix D; if the site contains a second-class particle we write a matrix A. 
The matrices D, E, A obey the algebraic rules 

DE - qED = (1 - q)(D + E) (10) 
DA-qAD = (1 -q)A (11) 
AE-qEA = (1 -q)A . (12) 

The only remaining condition to satisfy is that representations of this algebra may 
be found which give well-defined values for the traces appearing in (JTj). This may be 
achieved for q < 1 as follows. Let A be the diagonal, semi-infinite matrix 
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Then D,E may be chosen to be bidiagonal semi-infinite matrices 
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(14) 



(15) 



where 



Q 



n+l 



(16) 



We note that, as A has a geometric series for diagonal, the trace of A times a finite 
product of D and E matrices is not divergent. See [3] for further representations of 

(HMD). 

In principle, explicit formulae for the weights of each configuration may be obtained, 
either by using the algebraic rules or taking an advantage of an explicit representation 
of the matrices such as ( fT3HT5l) (see [20] or the review [3] for details of how these 
calculations are performed). 

In what follows, it will be useful to consider the matrices 5, e defined by 



S 



;i7) 

'18) 



D - 1 
e — E — 1 

which verify the algebraic relations 

Se - qe5 = (1 - q)t 
5 A = qA5 
Ae = qeA . 

The first equation (fl9|) is the (/-deformed harmonic oscillator algebra (see [15] for an 
introduction). Writing 5,e out explicitly we have 

\ 



(22) 
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In [28] the matrix product formulation of the stationary state of the iV-TASEP was 
presented, here we present the generalization to the iV-ASEP and give complete proofs 
that the solution indeed satisfies the stationarity condition. 

We first fix our notation. We denote the configuration of the system . . . , tl) 
by t. We use X^ to denote the matrix associated to the state t\ of site I in a system 
containing N species of particle. Thus, the subscript Ti = 0,1, ... N indicates to which 
species of particle the matrix is associated, and the superscript (N) indicates that 
there are N species of particles in the system. This is required because the matrix 
corresponding to a species will vary according to the total number of species in the 
system. The stationary probabilities become in the matrix product formulation 

■Tr[xWxW...XW] . (24) 



P( 



Z(Po,-..,^ 



AT; 



The matrices X^ are given by 



with 



X 



X, 



(N) 



(1) 



AT-1 

£• 

M=0 



L JM 



X 



(AT-1) 
M 



ioi < J < N 



(i) 



(25) 



(26) 



Note that we use a slightly different notation to [28] for the TASEP; with our notation, 
some of the will be equal to zero. 



The matrices aj^ are given by 



7 (A0 

l JM 

AN) 

l JM 

AN) 

l 0M 

AN) 

l JJ 

AN) 

l 00 

AN) 

l J0 

AN) 

l N0 





t ®{M-l) ^ e q t ®(N-M-l) 
A®(J-V g) 1®(N-J) 
1®(N-1) 

A® iJ -V (g, 5 i^^-^" 1 ) 

<N) 



for < M < J 
for < J < M < N 
for < M < N - 1 
for < J < N - 1 
for < J < M < N 
for < J < N - 1 



(27) 
(28) 
(29) 
(30) 
(31) 
(32) 
(33) 



Thus the matrices X K at level iV are composed of tensor products of dim^v = ( 2 ) 
fundamental matrices 5, e, A, or 1. For N — 1, dim^ = which implies the use of 
scalars for the single species case ff26l) (i.e. all configurations are equally likely). For 
N = 2, dim^r = 1, which implies the use of matrices as given by ffTBTfToT) . Let us now 
check that the 2-ASEP matrices are recovered. We find from equations ( 1271) . ( 1291) . ( 1301 ) 
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and (l33l) . that aig = 1, a™ = e, afrj = 5, a?) = 1, and a™ = A. Then, from the 



nx) ~~ "01 ~~ c ' "io ~~ ") "n ~~ ' ail> - 1 "20 ~~ ^- j-^h, 
recursion relation ([25]) we deduce that 



Xo ] = 4? + 4? = 1 + e = ^i (2) = 4o + a ii =t + S = D, (34) 
and XP = a§ = A , (35) 

as expected. 

(3) 

Let us now write out the case X = 3. Since dini3 = 3, the X K matrices are built 
with tensor products of three fundamental matrices. From the definitions (I27H33]) . we 
have a£J = 1 1, a$ = e 1, ajjj = 1 e, af } = 5 1, aSf = 1 1, af 2 } = 5 e, 
ago = A® 5, a 2 2 = ^4 1 and ag = A A Using the recursion relation fl25l) and the 
results X^ = £, x( 2) = D, xf ] = A, we obtain 

v(3) (3) „ v( 2 ) . (3) „ v(2) . (3) „ v(2) 

= a oo + a 01 ; Xj y + a 02 ' X^ 

= l0l0 J E + e0l0 J D + l0e0A (36) 

X{ 3) = ag X< 2) + ag> <g) x{ 2) + ag X® 

= <50l0£ + l0l0L> + 50e0,4 (37) 

Xf ] = ag + ag } Xf } = A 5 E + A 1 A (38) 

) = 4 3 } Xg 2) = A® A® E . (39) 

We notice in this example that the matrices xj^ have the same expression in terms of 
5, e and A as for TASEP [28J; the only difference lies in the deformation by q of the 
algebra (fT9H211 between 5, e and A. 

In later sections we shall establish the algebraic rules, satisfied by the matrices 
Xt^\ which generalize the quadratic algebra (fTUHTTl) of the N = 2 case. First we 



discuss how the recursive structure ( 1251) allows the stationary weights for an X-species 
system to be written in terms of those for an (X— l)-species system. 

4. Transfer matrix interpretation of the matrix product representation 



In this section, we show that the matrix product representation (1241) can be conveniently 
rewritten in terms of a transfer matrix acting on the configuration space. 

4-1. Definition of the Transfer Matrix 

We begin by noting that ( T2"4"|) and ( 1251) define the stationary weights for a system with 
X species of particle recursively in terms of the stationary weights for a system with 
X — 1 species: 

N-l 



Tr 



i E o ^[Hs«^r i, )-(<'£X"'- I) 

E^I(4S-42>K- 1) -^r 
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= E^"^]^"" 1 '"^] • (40) 

i 

We use the notation j = (ji, j L ) and i = (ii, where j is a configuration of a 

system with iV species of particle and i is a configuration of a system with N — 1 species 
of particle. Each ji can take values from to N and each ii can take values from to 
N — 1. The sum in ( l40i) is over all configurations i with iV — 1 species. 

Let us now introduce a notation for the configuration space which will be of use in 
the sequel. We call V^ N ' the N + 1 dimensional vector space corresponding to the N + 1 
possible states of a site for a system with N species. If we do not specify the number of 
particles of each species, the total configuration space of the system is v[ N) = (V^)® L . 
We denote a vector in this space corresponding to configuration j = (ji, j'l) by \j). 
The set of all (N + 1) L possible configuration vectors \ j) form a basis of V^. We denote 
the steady state eigenvector for a system containing N species of particle by \N) where 

\N) = J2 W (M- (41) 

j 

The stationary state weights W(j) (jHD are then given by 

W(j) = (j\N) . (42) 
Using this notation we can write relation (1401) in the form 

(J\N} = J2(M N) \1)(AN-1), (43) 

i 

where is the transfer matrix for a system with N species. The matrix element 

(j\T^ N '\i) can be thought of as representing a transition from configuration i in V^ N ~ lS> 
to configuration j in V^ N \ (This transition is, of course, not the same as a dynamical 
transition given by the Markov matrix.) The transfer matrix is used to express the 
stationary weights for a system with iV species linearly in terms of the weights for a 
system with N — 1 species. We identify from ( 140]) the elements of the transfer matrix 



Tf ) as 



(M N ^^Tr[a^>...^> L \ . (44) 

We can write the relation (143j) more simply as 

|iV) = Tf ) |iV-l), (45) 

and iterating we obtain 

|iV)=Tf ) ...T| 2) |l), (46) 

where the eigenvector of the system with only one species |1) is such that each 
configuration has the same weight. 

Let us now formalise the mathematical structure of T^ N \ The transfer matrix t£ 
is a (N + 1) L x N L rectangular matrix which takes a vector in V^ N ~^ and sends it 
to in V^ N \ It is expressed as a trace of a product of the local tensors al*\ From 



AN) (N) 
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are themselves tensorial products of elements of 
{5, e, A, 1}. Thus, the building blocks of the transfer matrix and indeed 



( I27H33I) . we observe that these a 
the set T 



the matrices Xj N>> are the four infinite matrices of the set T . We will call A the 
infinite dimensional space on which the elements of T act. The matrices af^ act on the 
auxiliary space A^ N > = A®^ -1 ^. They can be seen as the element JM (0 < J < A and 
< M < A — 1) of a rectangular (N + 1) x iV matrix a^ N \ For example, 
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5 <g> e 
A® 1 
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(47) 



4-2. Interpretation of the matrices 



In this subsection, we give an interpretation of the transfer matrices in terms 

of forbidden and allowed transitions between configurations of particles. Forbidden 
transitions correspond to matrix elements of (j\T^\i) that are equal to zero. We prove 
in Appendix B that the only nonzero matrix elements (j\T^\i) between an initial 
configuration i of particles of species between and N — 1 and a final configuration j 
of particles of species between and iV are characterized by the following rules: 

• at each site, a hole in the initial configuration i can remain a hole or become a 
particle of any class between 1 and A in the final configuration j. 

• at each site, a particle of class % between 1 and iV — 1 can become either a hole or 
a particle of class j between 1 and i (the class can only decrease). 

• there is global conservation of the number of particles of each class between 1 and 
N — 1. 



These rules are proven in Appendix B In particular, we note that the number of holes 
decreases (or stays the same) whereas the number of particles of class can only 
increase from none. Thus, the only way to create a particle of class N is to remove a 
hole and create a class A particle instead. The allowed local transitions for A = 2 and 
A = 3 along with the corresponding local tensors Ojf appearing in the transfer matrix 
element are illustrated in figures [1] and [2j 

To illustrate the utility of the transfer matrix, we work out some simple examples. 
First we consider the configuration (2, 1, 0) for L = 3 and A = 2. According to the rules 
for nonzero elements of the transfer matrix in the beginning of this subsection: the 2 at 
site 1 must have come from a 0; the 1 at site 2 could have come from a 1 or 0, the at 
site 3 could have come from a 1 or 0. Then the global constraint of a conserved number 
of 1 implies that the only 1-species configurations which have transitions to (2, 1,0) are 
(0, 1, 0) and (0, 0, 1). Using Figure [T]to construct the transfer matrix elements from the 
local tensors a,jf \ we find 



W(2,l,0) 



W(0,1, 0)Tr 
Tr [A] + Tr [A8e 



( 2 U 2 U 2 ) 

11 u 00 



a 20 a 



+ W(0,0,l)Tr 



(2) (2) (2) 
a 20 a l0 a 01 



(48) 
(49) 
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(2) 
a ji 



.4 



Figure 1. Allowed local transitions in the transfer matrix T L between a particle of 
class i (lower row, i = or i = 1) and a particle of class j (upper row, j = 0, j = 1 
or j = 2) at the same site. A particle of class is represented by a vertical dashed 
line, a particle of class 1 by a full line, and a particle of class 2 by a double line. The 

(2) 

corresponding local tensor appearing in the transfer matrix element is also shown. 



7 (3) 



5(8) 1 



A®5 



A® A 



A®t 



Figure 2. Allowed local transitions in the transfer matrix T£ between a particle of 
class i (lower row, i = 0, i = 1 or i = 2) and a particle of class j (upper row, j = 0, 
j = 1, j = 2 or j = 3) at the same site. A particle of class is represented by a vertical 
dashed line, a particle of class 1 by a full line, a particle of class 2 by a double line, 

(3) 

and a particle of class 3 by a triple line. The corresponding local tensor appearing 
in the transfer matrix element is also shown. 

where we have set 1^(0, 1,0) = W(0, 0, 1) = 1 (uniform measure for N — 1). This is to 
be compared with the known expression given by the matrices (Q 

W(2, 1, 0) = Tr [ADE] = Tr [A(l + 5)(1 + e)] =Tt[A + ASe] (50) 

where we have used the property that 5 and e must appear in equal numbers to have a 
non zero trace. 
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As another example we consider the configuration (3, 2, 1, 0) for L = 4 and N = 3. 
According to the rules for nonzero elements of the transfer matrix: the 3 at site 1 must 
have come from a 0; the 2 at site 2 could have come from a 2 or 0, the 1 at site 3 could 
have come from a 2,1 or and the at site 4 could have come from a 2,1 or 0. Then 
the global constraints of conserved numbers of 1 and 2 implies that the only 2-species 
configurations which have transitions to (3, 2, 1, 0) are (0, 2, 1, 0), (0, 2, 0, 1), (0, 0, 2, 1), 
(0, 0, 1, 2). Using Figure [2] to construct the transfer matrix elements yields 

1^(3,2,1,0)= (51) 
W(0, 2, 1, 0) Tr [afja^afla®] + W(0, 2, 0, 1) Tr 

+ W(0, 0, 2, 1) Tr \a®a®a®a$] + W(0, 0, 1, 2) Tr 



u 30 u 22 "10 "01 



(3) (3) (3) (3) 
°30 a 22 a ll a 02 



= W(0, 2, 1, 0) Tr [A 2 ] Tr [A] + W(0, 2, 0, 1) Tr [A 2 5e] Tr [A] 

+ W(0, 0, 2, 1) Tr [A 2 5e] Tr [A8e] + W(0, 0, 1, 2) Tr [A 2 ] Tr [A5e] . 

We used the fact that the trace of a tensorial product is equal to the product of the 
traces. 

4-3. Finiteness of the Trace Operation 

We have so far avoided the important question of the finiteness of the matrix product 
representation ( J24|) . Using previous results of this section and the appendices, we now 
prove that if there is at least one particle of each species, the expression (1241) for the 
stationary measure of the N species ASEP is finite. (In the case where there are zero 
particles of some species, this species can be removed from the problem by studying 
the corresponding system with N — 1 species instead of N.) We will focus on the case 
N = 3. For any configuration k with 3 species, the stationary weight is given by 

W(k) = (k\3) = ^^(fc|T| 3) |J)(J|Tf , (52) 

1 i 

where i is a configuration with 1 species and j a configuration with 2 species. Let us 
assume that there is at least one particle of each class in k. Then, as conserves the 
number of particles of class 2, the configurations j that give a nonzero contribution to 
(k\3) are such that j has also at least one particle of class 2. From the characterization 
of the nonzero elements of the transfer matrix, we must have that both i\ = 0, ji = 2, 
and ji' = 0, ky = 3 at some sites / and V between 1 and L so that i and j give a nonzero 
contribution to (k\3). We see, using the expression of the transfer matrix elements 



and the form ( 1331) of for N = 2 and = 3, that there will be at least one A in each 
trace of (j\T^\i) and (k\T^ \j) contributing to (fc|3). But one can calculate explicitly 
traces of products of elements of JF when there is at least one A in the product, and 



these traces are finite (see Appendix A ). This proves that (fc|3) is finite. The extension 
to arbitrary iV of this proof of the finiteness of the matrix product representation (1241) 
is straightforward. 
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4-4- Relation with Ferrari- Martin's construction for TASEP 

For the totally asymmetric case (q = 0), the transfer matrices encode Ferrari- 
Martin's multiline construction of the stationary weights j27j. We will focus on the 

(2) 

2-species stationary eigenstate constructed by T£ . We will see that the set of pairs of 
configurations i and j for which (j\T^\i) ^ is smaller than in the case q ^ 0. In order 
to know for which pairs of configurations i and j the matrix element [j\T^\i) is equal 
to zero, we refer to representation ( 1B.4j) of Appendix B which expresses an element of 
the transfer matrix as a product of iV — 1 traces of products of L fundamental matrices 
{5, e, A, 1}. Therefore, we need to know which products of elements of {5, e, A, 1} have 
a trace equal to zero. For the case of the TASEP, i.e. q — 0, the algebra ( fl9l) . ( 1201 . (12T1) 
between 5, e and A reduces to 

Se = 1 (53) 

cL4 = (54) 

Ae = . (55) 

One can see that for any product of / elements from {5, e}, w = w\...Wi, Aw A ^ if 
and only if both of the following conditions are true: 

• there are as many 5 and e in w. 

• for all m between 1 and there is at least as many e as 5 in w m ...wi. 

In that case, Aw A = AA = A. For example, AdSeedeA = A but AdeeddeA = because 
of an excess of 5 to the right. Then, each of the traces of products of elements from 
{5, e, A, 1} appearing in the transfer matrix element (1B.4j) will either be or 1 and 
consequently each transfer matrix element is either or 1. Then the transfer matrix 
relation fT43|) becomes an expression for the weight of a configuration in a system of 
species as a sum of weights for 'ancestor' configurations of a system of A^ — 1 species. 
The rules for selecting these ancestor configurations are precisely those given by Ferrari 
and Martin. 

We illustrate the equivalence for q = of the transfer matrix relation (TJ3]) with the 
Ferrari and Martin algorithm in the case N = 2, studied by Angel [29]. In this case 
(see Figure d]) an e in the transfer matrix element corresponds to a particle of class 1 
changing into a hole, a 5 corresponds to a hole changing into a particle of class 1 and 
an A to a hole changing into a particle of class 2. Therefore, for q = 0, (j\T^\i) ^ if 
and only if the configurations % and j are such that one can go from i to j by creating 
particles of class 2 at some of the unoccupied sites and by moving particles of class 1 
only to the left, forbidding them to cross class 2 particles: this is precisely the pushing 
procedure of Angel for 2 classes of particles which is a particular (2-line) case of the 
Ferrari-Martin AMine algorithm (see [28]). 

As a simple example of the distinction between the q = and q ^ cases let us 
consider the configuration (0,1,2) for L = 3 and N = 2. Constructing the transfer 
matrix elements as before, we find 

W(0, 1, 2) = W{0, 1, 0) Tr [agaMl + W{1, 0, 0) Tr fa^agag 
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= Tr [A] + Tr [e6A] , (56) 

where we have set W(0, 1,0) = W(l, 0, 0). In the case q = 0, Tr [eSA] = and there is 
only one contribution to W(0, 1,2). This concurs with the pushing procedure for this 
example which results in just one ancestor configuration (0,1,0). However, for g ^ 0, 



the second term in ( J56l) does contribute and using the result of Appendix A one finds 



Tr [eSA] 



l + q 



Tr [A] 



A further important observation for the 2-TASEP is that the matrix elements 
{j\T^\i) "decouple on the particles of class two", that is, each element of the transfer 

(2) 

matrix factorizes as a product over all the pairs of two consecutive particles of class 
2 in j of terms depending only on the holes and first class particles in i and j between 
the particles of class 2. This is because for q = the matrix A is a projector ffT3]) . 
However, this factorization property does not hold anymore for the general case q ^ 
as (Tl~3|) is no longer a projector. This fact makes it more difficult to find a combinatorial 
interpretation of the transfer matrix such as Ferrari-Martin multiline construction since 
the hops of the particles are less restricted than for TASEP. 



5. Quadratic algebra for the 



From the algebra for 5, e and A ([191421]) and the explicit form of the 

ajM (EMI, many 

quadratic relations for the be deduced. In the proof of the matrix product 

representation, to be presented in section El we will need two kinds of these relations: 
symmetries of a quadratic function of the a j^ under the exchange of the second indices 
of the a^l and relations allowing us to commute two cr 



"JM 



JM- 



5.1. Symmetry relations 

We have the following relation between commutators 



a 



(N) (N) 



JMi a KM' 



a 



(N) (N) 



JM' 



. a 



KM 



for J ^ and K ^ 0, 



(57) 



which indicates that 



a (N) a (N) 

U JM ' U KM' 



is symmetric under the exchange of M and M' for all 



M and M' . For some values of J, K, M and M', we have the stronger properties 



a {N) a (N) 



AN) (N) 



l JM 



,a^M for all J, K, and M, M' G {0} U [max( J + 1,K+ 1),N — 1] (58) 

(59) 



a jVi a jM' = a jM' a jIf f° r a ^ Ji M an d M' , 
which are also symmetries under the exchange between M and M' . 



5.2. Commutation relations 

In the following, we will also need to exchange a[^ M and a^2f /- ^ J = K, the symmetry 
relation (1591) can also be seen as a commutation relation for all value of J, M and M'. 
For J K, the exact form of the commutation relation will depend on M and M'. In 
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the case < J < K we partition the ensemble of couples (M, M' ) for which a j^ and 
a { p M , are defined and non-zero (that is M e {0} ll[J,N— 1] and M' e {0} U [if, iV - 1]) 
into 12 sectors as 

M = J < M < K M = K 

1 4 7 

2 5 8 

3 6 9 



W = 
M' — K 
K < W < N - 1 



K < M < N — 1 
10 
11 
12 



.(60) 



Then, we have the following commutation relations between a^j and a^^, for < J < 
K: 



a (N) a {N) 

U JM U KM' 

a (N) a iN) 

U JM U KM' 

a (N) a (N) 

JM KM' 



(N) (TV) 
1 a KM' a JM 

a {N) a {N) 

U KM' U JM 



in sectors 1, 2, 3, 8, 10, 11 and 12 



in sectors 4, 5, 6 



a 



(N) (N) 
KM' U JM 



q)a 



KM JM' 



in sectors 7 and 9 . 



(61) 
(62) 
(63) 



In the case K = < J, the commutation relations take a similar form (IC.2HC.4]) 



6. Proof of the matrix product representation 

In this section, we prove that the matrix product expression (J24l) gives the stationary 
state eigenvector of the Markov matrix. 

6. 1 . l Hat ' matrices 



The Markov matrix of the system with N classes of particles can be written in terms 
of the local (N + l) 2 by (N + l) 2 matrices M$ +1 which encode the rates at which the 
particles hop between site k and k + 1 . 

L 

M = J2 ® M l N k+i ® t^ L ~ k -^ . (64) 

k=l 

For a model for which the rates do not depend on the site such as the one we are 
discussing, the local matrices do not depend on the site: M^ +l = M^J . For example, 
for N = 2, in the basis (11,12,10,21,22,20,01,02,00) 

/ 

-1 • q 



M, 



(2) 



loc 



-Q 



\ 



(65) 
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The dots represent matrix elements equal to zero. The matrix product expression 
gives the stationary state eigenvector if for all configuration j with a particle of class jk 
at site k 



x {N) x (N) 



0. 



(66) 



or in terms of the local jump matrix 

L TV 



E E 



M 



(TV) 
loc 



Tr 



fc=l ifc.i 



Jkjk+l, l k l k+l 



Y (N) Y (N) Y (N) Y (N) Y (N) 



.X 



(TV) 



. (67) 



This equation will be satisfied if there exists some additional 'hat' matrices [30], [31] 



X, 



(TV) 



xff such that 



TV 

E 

i,i'=0 



M, 



(TV) 



ioc 



^(TV)^(TV) 



(6* 



leading to a cancellation of all the terms two by two. Knowing the form of the local 
matrix , the previous equation ( 1681) can be rewritten as: 

for < J < N 

, 



X (N) £(N) 



for 0<J<K<N or 



K < J <N 



(69) 



x^x^ 



AN) 
L J 



^(TV)^(TV) 



(70) 
(71) 



Again, the fact that the case = K < J < N is singled out comes from our choice to 
give the holes the index instead of N + 1 as the particle hierarchy would have required. 
Subtracting ( 1701) from ( JTTi) we also have the additional relation valid for all values of J 
and K, 



y(N) Y (N) 



(TV) V>(TV) 



X£>,X 



. 



(72) 



This equation holds even for J = K because of equation (1691) . It tells us that 
Xj N \ must be antisymmetric under the exchange between J and K. 

In the next subsection we will prove by induction on iV that the following X^ N > 
matrices and the X^ N ' matrices defined in equation (125|) verify equations (169]) and (170]) : 



X, 



X 



(TV) 



(TV) 



-(l-q)X l 



(TV) 



TV-1 

E' 

M=0 



AN) 
l 0M 



X 



(TV-1) 
M 



TV-1 

E' 

M=0 



VM 



with Xq 1 = 1 and X\ l> = q. The a V j^ are still given by equations (127])- (133]). In contrast 
with the Xj N ^ matrices (125]) . the expression for the X 
exactly the same as for TASEP: q enters both ( 1731) and the value of X{ 



^(TV-l) 



AN) 



for 1 < J < tV, 



(73) 



(74) 



(TV) 
J 



in terms of the is not 
(i) 
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We first check the simple cases N = 1 and N = 2. Using X^ = 1, X} 



X^ = 1 and X\ L> = q, we see that relations ( |69i) . (1701) and (|7T|l are trivially satisfied for 



N = 1. For X = 2 we find 



(2) 



"(2) 



,(2) 



X 



(1) 



,( 2 ) 



'01 



— (1 — g)X -r u, 00 

- (1 - q)E + ql + e = qE - (1 - q)t 



(i) 



(2) 



t (2) 

no 



X 



(i) 



(2) 



g5 + 1 = qD + (1 



if 



X® = < ; ® X£> = qA 



,(2) 



"(1) 



(2) 

Clearly, because XA 



E, X[ 2) = D and X^ 2) 



A, the commutation relations (1691) are 
satisfied. One can check, for example, that 

Xf ^Xf > - Xf 'I? = (QD + (1 - - D(qE - (1 - g)l) 

= (1 -<?)(£> + £) 

W " X { N) Xr = E(qD + (1 - q)t) - (qE - (1 - 

= (1 + 

so that (J7U1) and (1711) are satisfied for J = 1, X = and the other cases are similarly 
verified. 

In the next subsection, we prove the three quadratic relations f l69|) . (1701) and (1711) 
which are sufficient for the matrix product expression ( 1241) to give the stationary state 
eigenvector of the Markov matrix. For each relation, we must distinguish the case where 



some indices are equal to zero (which will be proved in Appendix C ), giving six relations 
to prove. The proofs all have a similar pattern. They are done by induction on N and 
rely heavily on the symmetry and commutation relations of section [51 



6.2. Proof of equation ( fgPjj 
If J 7^ we have 



N-l N-l 



(75) 



M=0 M'=0 



I i a JM a JM' 



(A r_l) A(JV-l) 
M M' 



AN) (N) 
l JM' u JM 



M' M 



Using the quadratic relation (l59l) to exchange fljjj/ and aj A j, we get 

N-l N-l 



X (N) X (N) 



M=0 M'=0 



M > M' 



(76) 



Using the fact that o!jM a jM' * s symmetric in M and M' (159|) . and that from 

is antisymmetric in M and M' by induction, we find that 



Y (N-1) y(N-l) 



X (N) X (N) 



= 0. 
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For J = 0, we have an extra term compared to (J75J): 

[xr.-d- q )xr\ 



A ! A 



(77) 



N-l N-l 



1^ 1^ [[ a 0M a 0M') ® [ A M A M' 



A _ (a (N) a {N) 

\ u QM' u 0M 



MN-1) Y (N-1) 
M' M 



M=0 M'=0 



But this extra term, being the commutator of Xq N ^ with itself, vanishes and (1771) reduces 
to (175]) in the case J = giving also 



A o > A o 



0. 



In the rest of the proofs, we will use the following convention to lighten the notation: 
Xm and Xm ^ will be written respectively Xm and Xm, while the a/j M will just be 
written o/m- 

6.3. Proof of equation (7§ (0 < J < K) 
Let us define 



A = X^ N) X^ ] - qX ( ^X { p - X^X^ + X { j N) X { p 



(71 



We want to show that A = 0. We have: 



A = 



^ ^ [{ajMdKM') ® (X M X M >) - q(a K M>aj M ) ® (X M /X M ) (79) 

j\/g{o}u[j,tv-i] M'e{o}u[i<r,JV-i] 

- {ajuCLKM') ® (^j/4f) + (a JM a KM >) <g> (X M X M ,)] . 

We will cut the double sum into 4 parts and write A — A1 + A2 + A3 + A4, gathering 
sectors from the partition (160]) . ^ will be made of the sectors 4, 5 and 6 of (160|) . ^4 2 of 
the sectors 1, 3, 10 and 12, A3 of the sectors 2, 7, 9 and 11, and A4 of the sector 8: 

Ai : M G [J, K - 1] and M' G {0} U [A", iV - 1] (80) 
A 2 : M G {0} U [K + 1, AT- 1] and M' G {0} U [K + 1,N- 1] (81) 

M G {0} U [i<f + 1, iV — 1] and M' = K 
3 '' M = K and M' G {0} U [A + 1,7V- 1] 

A : M = K and M' = K . 
We will now show that ,4i, A, ^.3 and A4 are all equal to zero. 



(82) 
(83) 



We begin with Ai and use the commutation relation (I6"2"]) . We get 



K-l 



A\ — {ojjmO'KM')® \ XmXm' — qXu'Xu — XmXm* + ^m^m'| -(84) 

M=J M'e{0}U[A _ ,AT-l] 

By induction, XmXm' — qX M iX M — X M X M i + X M X M i vanishes (1701) . and thus A\ = 0. 



For ^2, using the commutation relation (16T]) . we obtain 

v4 2 = ^ (ajua-KM') ® \ [Xm,Xm'] — X M X M i + X M X M i 

M,M>£{0}U[K+1,N-1] 



55) 
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From (I5S|) . cbjM a KM' is symmetric under the exchange of M and M' while [X M ,X M >] is 
antisymmetric, as well as — XmXm* + XmXm< which follows from a rewriting of f!72|) . 
by induction. This gives .4.2 = 0. 

For A3, using the commutation relations ( 16TT) and ( |63l) . we have 

A 3 = ^ [{o-jkO-km 1 ) ® \ X K X M i — qX M iX K — X K X M i + X K X M /\ 

M'e{0}u[K+l,N-l] 

- q(l - q) (a KK a JM i) ® (X M >X K )] 

+ (O'JMO'Kk) ® ^M^/f — -XaT-^m — XmXk + A'mA'k) (86) 

Me{o}u[ftr+i,jv-i] 

The first sum corresponds to setting M = K and the second to M' = By induction, 
XkXm' — qXM'Xx — XkXm' + XkXm 1 = in the first sum ( TfOj) . With the help of the 
commutation relation ( IBTl) in sectors 2 and 11 for clkkQ'JM' , we obtain, after renaming 
the dummy variable M' to M 

A3 = (cijmcikk) ® ^^Cm^Cr: — X K X M — X M X K + j . (87) 

Me{o}u[Jf+i,JV-i] 

By induction, qXuXx — XkXm — XmXk + XmXk = from (JTB . as only M = and 
M > K contribute to the sum. A3 is therefore equal to 0. 

Finally, for Aa, we have 

Aa = (a JK a KK - qa K Ka JK ) <8> (X K X K ) + (cijk&kk) <%> (-X K X K + X K X K J .(88) 

In the first term cljko^kk — qdKKdjK = because of the commutation relation ( I6"T1) . 
and in the last term — X K X K + X K X K = by induction because of fl69|) . Thus *A 4 = 
which proves (1701) in the case < J < K. 

6.4. Proof of equation (7$ (0 < J < K) 
Let us define 

B = X™X™ - qX\Px^ - xP&f> + jt™X™ (89) 
We want to show that B = 0. We have 

B = ^2 \(flJM a KM ,s ) ® P^m^m') - q(a KM iaj M ) ® (X M /X M ) (90) 

A/e{0}U[J,V-l] M'6{0}U[Jf,Ar-l] 

— (clkm'Cijm) <S> (Xm'^m) + (aRM'ttjAi) <8> (X A //Xm)] ■ 

Again, we cut the double sum into 4 parts and write B = B± + B 2 + B 3 + £> 4 , gathering 
sectors from the partition ( JBDI) . We use the same sectors as in the previous proof ( j80T 
1831) . At this point, we have to calculate £>i, #2, S3 and £> 4 . Using the same arguments 
as in subsection 16.31 we find that they are all equal to zero, which proves (ITT]) in the 
case < J < K. 
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7. Conclusion 

A solution for the stationary state of the multispecies TASEP was first proposed by 
Ferrari and Martin in [27]. Their solution was expressed in terms of numbers of 
configurations of a multiline queueing process. A matrix product representation of this 
solution was then given in [28] for the case of periodic boundary conditions, making the 
link with several works on one dimensional exclusion processes in the physics literature. 
In this paper, we have extended this matrix product solution to the previously unsolved 
problem of the multispecies partially asymmetric exclusion process on a ring. In this 
case, there is no known analogue to Ferrari-Martin queueing construction. 

The mathematical structure of the matrix product solution for the stationary state 
reveals several interesting features. First, the matrices are defined in a recursive fashion 
(see eq. ( )25il ) using auxiliary matrices a-j , and are ultimately built as tensor products 
of the fundamental matrices {5, e, A, 1} used in the N = 2 solution. Second, the matrices 
obey quadratic relations involving additional hat matrices and these relations generalize 
the quadratic algebra of the N = 2 case (jTOHTTj) . Such relations have only been verified 
before in a few cases (see [3] section 9 for a discussion). In our solution the key to 
satisfying these conditions lies in the algebraic properties of the auxiliary matrices 
which we presented in section [5] 

The recursive structure of the solution allows us to construct a transfer matrix 
relating the stationary weight of a configuration with N species of particles to the 
weights of configurations with iV — 1 species. For the case q = the transfer matrix 
recovers the algorithm of Ferrari and Martin, whereas for q ^ the structure is more 
complicated. 

In this paper we have not attempted to calculate physical quantities of interest, 
such as correlation functions. Such calculations begin with the computation of the 
normalisation Z(P , . . . , Pn), defined in ( 12^1) . for all system sizes and particle numbers. 
Even for the two species case this computation is not easy and to our knowledge has only 
been carried out for the totally asymmetric case |20j. The transfer matrix construction 
of section |4] may provide a formalism for the computation of the normalisation for partial 
asymmetry and general N. 

There is a well known correspondence between the matrix product representation 
for the one species ASEP with open boundaries and the two species ASEP on a ring: the 
matrices corresponding to the holes and to the first class particles are the same in both 
cases, but the matrix corresponding to the second class particles on the ring becomes 
associated with the boundaries for the open system. This correspondence should be 
investigated further in the case of general iV species systems. Another interesting 
extension is the case of systems with different rates for the different classes of particles, 
which contains in particular the ABC model [32] . 

The multispecies asymmetric exclusion process shares with a small number of 
statistical physics models the significant property of being integrable. In particular, 
this means that its Markov matrix can be diagonalized using the (nested) Bethe 
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Ansatz. The matrix product representation of the stationary state and the Bethe Ansatz 
[331 El ESI ESI HOl EH M, U M, E] are two of the most used techniques used to obtain 
exact results for the ASEP. Thanks to its rich structure, the solution of the multispecies 
ASEP might help to understand more precisely their relationship. 
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Appendix A. Traces of products of S, e and A 

In this appendix, we show that the product of elements of the set T = {5, e, A, 1} with 
at least one A has a nonzero trace for q ^ if and only if there is the same number of 8 
and e in the product. (The case q = has already been dealt with in section FOl ) We 
use an explicit calculation of one of these traces. For g^Owe can first use the relations 
( 1201 I2TI) to rearrange any product of elements of T by putting all the A to the left. The 
result will be proportional to r Yx{A p w) where w is a product of elements from the set 
{5, e}. We call r the number of 5 in w and s the number of e. Using again (1201) and (12T1) 
to commute all the A matrices to the right, we get 



Using the cyclicity of the trace, we see that if p > 1 and r ^ s, Tr(A p w) must be equal 



If p > 1 but r = s, the algebra (Tl9l) between 5 and e allows one to rewrite wasa 
linear combination with nonnegative coefficients (for q between and 1) of terms of the 
form e r '5 r ' with r' < r. We want to show that Tr(A p e r ' 5 r ') > in order to prove that 
Tt(A p w) if p > 1 and r = s. 

Let us define 



We now show how fi can be calculated easily by recursion on r. Using first the 
deformed commutator ffT9]) between 5 and e once, we get 



Tt(A p w) = q p{s ~ r) Ti(wA p ) . 



(A.l) 



to 0. 



fl p) = Tr(A p 5 r e r ) . 



(A.2) 



= Tr(A p 5 r (qe5 + (1 - q))e r ) = (1 - g)/W + gTr(^W) . 



(A.3) 



Using ( fl9i) r more times to push the rightmost 5 to the right we obtain 

/ r W = ... = (1 - g)(l + g + ... + g fe )/ r (p) + g fc+1 Tr(yWe fc+1 <k r " fc ) 



= ... = (1 - q r+1 )f^ + <f +1 Tr(A p 5 r e r+1 5) . 



(A.4) 




(A.5) 



We have found the recurrence relation 



1 nT-TX 



(A.6) 
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which gives 



Tr A p 
r + p 
V 



J q 



Here 



is the q-deformed binomial coefficient (see e.g. |15j ) defined as 



[a] q [a- 



q->l 



([6] a [6-l] ff ...[l] 9 )([a-6] ff [o-6-l] 9 ...[l] ff ) 



the q-deformed numbers [a] 5 being defined by 
r i 1 " <f 



l + g + ... 



Also, using representation (fl3|) we may evaluate 

1 



<r 1 ^ a . 



1 - qP 



(A.7) 



(A.J 



(A.9) 



(A.10) 



Thus, the f~r are positive for all r' up to the factor TyA p , which means that 
Ty(A p w) cannot be equal to zero if p > 1 and r — s. We have found a necessary and 
sufficient condition for a product of elements from T to be different from 0. 



Appendix B. Proof of the characterization of the Transfer Matrix 

We now prove the characterization of the nonzero matrix elements of the transfer matrix 

(j\T^ N '\i) is a trace of products of . Because 



T^ N \ given in section I4~2l From (J4 



of (1271) . if there exists a site I between 1 and L such that < i\ < ji, then (j\T^\i) = 0. 
Thus, for all /, we must have either i[ = or ji < ii for the matrix element of 
to be nonzero. In terms of the particles in the configurations i and j, this means that 
a hole in i can become a particle of any class between and iV by application of the 
transfer matrix, and that a particle of class i\ > 1 can become a particle of class ji only 
if < ji < i;. For example, for N = 2 (fig. [T|), only the five following transitions are 
allowed in the transfer matrix : 0— »0, 0— »1,0— >2, 1— »0 and 1 — > 1. The transition 
1 — ^ 2 is forbidden. For N = 3 (fig. [2]), 9 transitions are allowed: 0— >(), 0— >1,0— >2, 
-»• 3, 1 -»• 0, 1 -> 1, 2 -»■ 0, 2 -»■ 1, 2 -»■ 2. The transitions 1 -> 2, 1 -> 3 and 2 -> 3 
are forbidden. 

— * — * 

But this local constraint on the classes of the particles in i and j at each site 
/ does not characterize completely the nonzero matrix elements of we will see 

that there is also a non-local constraint on % and j. We observe that the expression 



for the transfer matrix T- 



(TV) 



can be simplified further by noting that the are 



themselves tensorial products of iV — 1 elements of T 
the notation 



{6, e, A, 1} (J2ZH335. Introducing 



(B.i; 
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we have 



(JV.l) (iV,m ( (N,N-1) (N,N-1)\ mQ v 

the a^'^ being elements of JF. Thus, we have for the transfer matrix element 

<J|7f>|!> = IV (a£"... S;») x ... x TY (af^.^) (B.3) 

jV-1 

(B.4) 



niv 



n 

.1=1 



Representation flB.41) shows that an element of the transfer matrix can be written as 
a product of iV — 1 traces of products of L fundamental matrices {5, e, A, 1}. In the 



following, we shall call Tr Ylt=i a ^i^ the t race °f O'I^l V; |*)- As is expressed 
in terms of traces of elements of JF, we have to study which elements of T have zero 
trace and which have nonzero trace in order to determine which are the nonzero matrix 
elements of T^ N \ In Appendix A we showed, using the quadratic algebra ( JT9142TI) . that 
the trace of a product of elements of T with at least one A is nonzero (for q ^ 0) if and 
only if the number of 5 is the same as the number of e in the product. In terms of the 
matrix elements of the transfer matrix , this means that each of the A — 1 traces 
of (j\T^ \i) must contain the same number of 5 and e as we now show. 

The condition of having at least one A in each trace is automatically verified if we 
choose the configuration j such that it contains at least one particle of class A" because 
a factor appears (133|) . From the expressions (I27H331) for the , we see that a 
5 appears in the rth trace of (j\T^\i) for all I such that ji = r and either %\ = or 
%x > r. This corresponds to sites at which a hole or a particle of class strictly larger 
than r is replaced by a particle of class r. On the other hand, an e appears in the rth 
trace of (j\T^ \t) in the expression flB.4j) for all I such that i\ = r and < ji < r. 
This corresponds to sites at which a particle of class r becomes a hole or a particle of 
class strictly lower than r. For example for A = 2, 5 corresponds to — ► 1 and e to 
1 — > (see fig. [1]). For N = 3, a 5 appears in the first trace for — ► 1 and 2 — > 1 and 
an e appears for 1 — > 0. And there is a 5 in the second trace for 0^2 and an e for 
2^0 and 2 — * 1 (see fig. [2]). In fig. IB 1 L we draw all the transitions involving a particle 
of class 2 when A = 4 (the ones which are forbidden are crossed out), along with the 
corresponding values for cif^ (second trace). To summarize, in the rth trace: 

• a 6 appears when a particle of class 0, r + 1, r + 2, ... , A — 1 in z is transformed 
in a particle of class r at the same site in j. 

• an e appears when a particle of class of class r in % is transformed in a particle of 
class 0, 1, 2, ... , r — 1 at the same site in j. 

Therefore, the requirement that the number of 5 and e is the same in the rth trace of 
in (IB. 41) implies that the number of particles of class r between 1 and A — 1 is 
conserved by the action of the transfer matrix T^ N \ But we emphasize that neither the 



< N )\7\ A o rp(N) 
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Figure Bl. All the situations involving a particle of class r = 2 for N = 4. A particle 
of class is represented by a vertical dashed line, a particle of class 1 by a full line, a 
particle of class 2 by a double line, and a particle of class 3 or 4 by a triple line. The 
crossed-out diagrams correspond to situations which are forbidden by ((27|) . 



number of holes nor the number of particles of class iV is conserved: the number of holes 
decreases while the number of particles of class N increases. This concludes the proof 
of the characterization of the nonzero matrix elements of the transfer matrix T^. 

Appendix C. Proofs of equations ( TTOl) and STB) for K = < J 

The commutation relations between and a^M' were given in section [5] for the case 
< J < K. They were used in section |6]to prove equations (TTDjl and (ITIj) for < J < K. 
In this appendix, we will write the commutation relations between and a^], and 
use them to prove equations (1701) and (JTT1) for K — < J. 

Appendix C.l. Commutation relations 

We partition the ensemble of couples (M, M') for which and a^, are defined and 
non-zero (that is M e {0} U[J,N- 1] and M' G [0, N - 1]) into 12 sectors as 

M' = < M' < J M' = J J < M' < N - 1 
M = T 4 7 10 
M = J 2 5 8 TT 
J<M<iV-l 3 6 9 12 

Then, we have the following commutation relations between and o!^},: 

a jM a oM> = Q a oM' a jM in sectors 4, 5, 6 (C.2) 

a jM a oM' = %M' a JM i R sectors 1, 2, 3, 8, 10, 11 and 12 (C.3) 

a jM a oS' = Q a oM> a j N M + C 1 - v) a 0M ai jM> in sectors 7 and 9 . (C.4) 



.(CI) 



Appendix C. 2. Proof of equation (T70| j (K = < J) 
Let us define 

7 = xf ] xi N) - qX^X^ - X^X^ + X?)XW . (C.5) 



A 



3 



(C.9) 
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We want to show that A = 0. With the extra term —(1 — q)XQ N ^ that comes in the 
recursive expression for X^ N \ the expression for A rewrites 

N-l 

A= ^2 ^2 ll( a JMaoM>) <8> (X M X MI ) - q(a 0M >aj M ) <S> (X M >X M ) (C.6) 

A/G{0}U[J,7V-1] M'=0 

- (a JM a 0M >) <g> (X M X M >) + (aj M a 0M >) ® {X M X M >)] . 

Like in the case < J < K, we will cut the double sum into 4 parts and write 
A = A\ + Ai + A3 + Ai, gathering sectors from the partition (IClft . A\ will be made 
of the sectors 4, 5 and 6 of fIC.lj) . Ai of the sectors 1, 3, 10 and 12, A3 of the sectors 2, 
7, 9 and 11, and A4 of the sector 8: 

Ai: M G {0}U [J,N- 1] and M' G [1, J- 1] (C.7) 

A 2 : M G {0} U [J + 1, TV — 1] and M' G {0} U [J + 1, TV — 1] (C.8) 

M = J and M' G {0} U [J+ 1, iV - 1] 
M G {0} U [J + 1, N - 1] and M' = J 

At: M=JandM' = J. (CIO) 

We will now show that A\, A2, A3 and A4 are all equal to zero. 

We begin with Ai and use the commutation relation ( 1C.2I) . We get 

JV-l J-l 

A\ = ^2 ( a J Ma oM') ® [qXmXm 1 — Xm'Xm — XmXm' + Xm-Xm'J ■ (C.ll) 

M=0 M'=l 

But, from (ITTj) . qX M X M i—X M iX M —X M X M i+X M X M i = by induction. Thus, A\ = 0. 

For A2, using the commutation relation (IC.3I) . we obtain 

A2 = ^2 (ajM a oM') ® [Xm,Xm'] — XmXm 1 + XmXm'^ ■ (C.12) 

M,M f e{0}u[J+l,N-l] 

From (|58|) . ajM&sM 1 is symmetric under the exchange of M and M' while LY M ,X A/ /] is 
antisymmetric, as well as —X M X M > + X M X M i because of ( 1721 by induction. This gives 

A 2 = 0. 



For A3, using the commutation relations (1C.3I) and (1C4j) . we have 



A 3 = (ajjaoM') <S> (q [Xj, X M ,} - XjX M , + XjX A 

M'e{o}u[J+i,iV-i] 

[( a JM«oj) ® (qX M Xj - XjX M - X M Xj + X M Xj\ 



Me{o}u[J+i,JV-i] 

+ (l-q) (a 0M ajj) <g> (XjX M )](C.13) 

By induction, qXuXj — XjXm — XmXj + XmXj = in the second sum ( ITT]) . With 
the help of the commutation relation (1C.3I) in sectors 2 and 11 for aQMQ-jj, we obtain, 
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after renaming the dummy variable M' to M 

A 3 = Yl (ajjaoAi) ® (XjX M - qX M Xj - XjX M + XjXm) . (C.14) 

A/e{0}U[J+l,7V-l] 

From (1701) . XjXm — oXuXj — XjXm + XjXm = by induction. A3 is then equal to 
0. 

Finally, for A4, 

A~4 = q (ajja 0J - a 0J ajj) ® (XjXj) + (ajja j) ® (-XjXj + XjXj) . (C.15) 

In the first term ajjaoj — clqjcljj = because of the commutation relation (lC.3j) . and 
in the last term —XjXj + XjXj = by induction (16"9"1) . Thus At = 0, which proves 
(170]) in the case X = < J. 

Appendix C.3. Proof of equation (71\ ) (K — < J) 
Let us define 

B = X^X^ _ qX^xW - XW*™ + X^X^ . (C.16) 

We want to show that B = 0. With the extra term — (1 — q)XQ N ^ that comes in the 
recursive expression for X ( N \ the expression for B rewrites 

N-l 

B= 22 ^ [(ajMOoM') ® (X M X M >) - (aoM'ajM) ® {X M ,X M ) (C17) 

M6{0}U[J,JV-1] M'=0 

- (a M>aj M ) <S> (X M 'X M ) + (oom'^jm) <8> (X M /X A/ )] . 

Again, we cut the double sum into 4 parts and write B = B\ + B 2 + B 3 + £> 4 , gathering 
sectors from the partition ( IC.lj) . B\ will be made of the sectors 4, 5 and 6 of ( 1C.1I) . B2 
of the sectors 1, 3, 10 and 12, B$ of the sectors 2, 7, 9 and 11, and B4 of the sector 8. 



Using the same arguments as in Appendix C.2 we can show that B\, B2, B3 and £>4 are 



all equal to zero, which proves (ITTj) in the case K = < J 
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